Separation of plastic deformations in polymers 
based on elements of general nonlinear theory 
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We report a method for describing plasticity in a broad class of amorphous materials. The method is based on 
nonlinear (geometric) deformation theory allowing the separation of the plastic deformation from the general 
deformation tensor. This separation allows an adequate pattern of thermodynamical phenomena for plastic 
deformations in polymers to be constructed. A parameter, 8 describing the stress relaxation rate of the material 
is introduced within the frame of this approach. Additionally, several experimental configurations to measure 
this parameter are discussed. 
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Plastic and elastic deformation plays a critical role in the 
description of physical phenomena in Earthis crustal defor- 
mations d, transfer of slurry, powdery and granular materi- 
als through pipelines |2j, and, especially, in polymer dynam- 
ics J2I. Although the theory of elastic deformation including 
elastic equilibrium, deflection and torsion of rods, bending of 
plates and shells is a classical and well developed field, a com- 
pleted theory applied to plastic media is yet to be developed. 

The theories describing microscopic plastic deformations 
in metals is based on a born-death dislocation mechanism pro- 
ducing a shift of crystallographic grid. This have been devel- 
oped by E. Orowan, M. Polanyi, G. I. Taylor, J. M. Burg- 
ers, F. C. Frank, W. T. Read, R. E. Pierels, P. B. Hirsch, 
W. C. Dash, Yu. A. Osipyan et al (0j). In solids, the shear 
transformation zone theory (STZ) has been developed by Falk 
and Langer 0|. However, a complete theory of plastic 
deformation in polymer materials has not developed so far. 
Physical behavior of polymers at the nanoscale is important 
from the fundamental point of view to understand deforma- 
tion processes, evolution of homogeneous or heterogeneous 
nanostructures, and the thermal history of monomer architec- 
ture Qlal. Understanding the deformation behavior in poly- 
mers would eventiually produce novel nanostructure forma- 
tion techniques based on nano-deformations. 

Polymers are probably most suitable materials in which 
plastic deformation can be studied experimentally. Recently, 
Lyuksyutov and Vaia with co-authors reported nanopattering 
technique based on localized Joule heating of a thin polymer 
films |9j 1 1 Oil . A biased atomic force microscope ( AFM) tip 
produces an electric current flow through the polymer film 
resulting in localized Joule heating of the polymer above 
its glass transition temperature due to electronic breakdown 
through the film. Polarization and electrostatic attraction of 
softened polymer toward the AFM tip in the presence of a 
strong (10 9 -10 10 Vm _1 ) non-uniform electric field produce 
raised or depressed nanostructures (10-50 nm width, and 0.1- 
100 nm height) in a broad class of polymers of different 



physical-chemical properties. This technique, named AFM- 
based electrostatic nanolithography (AFMEN), can be applied 
to the study of plastic deformations since both a liquid and a 
plastic solid polymer phases coexist during the process. The 
breakdown during AFMEN is a critical factor causing film 
softening, and polymer mass transfer as a result. Recent ex- 
perimental data indeed produce the evidence of nanostructure 
formation in polymer that cannot be explained by the elec- 
tronic breakdown. The most likely reason for polymer nanos- 
tructure formation when an AFM tip of 20-50 nm in diameter 
moves above the surface is the plastic deformation of polymer 
molecules through triboelectrification mechanism. No elec- 
tronic breakdown is required to deform the polymer surface in 
this case. An exact analytical solution, based on the method 
of images, has been obtained for the description of the elec- 
tric field between an atomic force microscope (AFM) tip and 
a thin dielectric polymer film (30 nm thick) spin-coated on a 
conductive substrate. Three different tip shapes are found to 
produce electrostatic pressure above the plasticity threshold 
in the polymers up to 50 MPa Hill . Should such a technique, 
based on plastic deformations in polymer materials, be devel- 
oped further, it would create an alternative to existing nanopat- 
tering techniques effective tool, for patterning on nanoscale. 
There has been a lot of activity in in last few years in devel- 
oping these techniques including: 1) direct resist lithography 
developed by Schaeffer et al II 211 based on the competition of 
Van der Waal's and Laplace forces on polymer-air interface in 
strong electric field; and 2) hierarchic nanostructure formation 
based on electrodynamic instability in bilayer polymer films 
developed by Russell et al JT3I1 . The only industrial proto- 
type for nanostructure formation in polymers called MILLI- 
PEDE developed by Vettinger et al [ 14] is based on thermal- 
mechanical lithography developed by Mamin and Rugar year- 
lier in the 90s. fij- The authors of the MILLIPEDE project 
predict a replacement of ferromagnetic memories with based 
on polymers in the next 20 years. 

Experimental verification of the existence of plastic defor- 
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mation on nanoscale in AFM-tip-polymer-metal system lHUl 
requires a theoretical model to describe them explicitly. An 
important question is how to describe a polymer surface un- 
dergoing deformations near the glass transition point when 
two phases exist. There are two approaches for describ- 
ing polymers in the liquid-solid phase deformed under ex- 
ternal forces. The first, would be based on intensive math- 
ematical description of this process through solution of the 
Navier-Stokes equation for a steady flow of non-Newtonian 
incompressible liquid with non-slip boundary conditions flfUl . 
Although this approach may be ultimately correct, it lacks 
generality for the description of polymer deformation at the 
nanoscale. The second approach would be based on time- 
dependent geometry of elastic and plastic deformations using 
the elements of differential geometry and tensor analysis. In 
this paper we consider this approach. 

The goal of this letter is the separation of plastic deforma- 
tions within a general nonlinear deformation tensor. A general 
method for the separation would be to build into a framework 
of balance equations traditionally used in the description of 
dynamics and thermodynamics of moving continuous media. 
The model we present, would be useful for the description of 
solid-liquid substances including organic molecules like DNA 
and long polymer molecules. 

The technique is based on the tensor calculus associated 
with curvilinear coordinates system. The basic concept for 
moving continuous media description is presented through the 
deformation map. The map transforms the point (y 1 , y 2 , y 3 ) 
of the non-deformed medium to its current actual position 
(y\ y 2 , y 3 ). Direct and inverse deformation maps can be pre- 
sented through the following sets of functions: 



y l = y l {y\f,f), f y 1 = y\y\ y\ y% 
v 2 =y 2 (y\f,f), \f = f{y 1 ,y 2 ,y% d; 

y 3 = y 3 (y\j/ 2 ,y 3 ), 



^2 _ r ;2, ,,! „2 „:i , 

f = f(y\y 2 ,y 3 ), 



The direct and inverse transformations of the deformation map 
are presented in Figure 1 [2 The partial derivatives of the map- 
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FIG. 1: Schematic presentation of the direct and inverse transforma- 
tions on the deformation map. A point is transformed 

ping functions Q define two Jacobi matrices 5* = dy i jdyi 
and T l : — dy 1 jdyK The non-linear deformation tensor can be 



defined using one of them as pesented below: 

3 3 

G l3 =Y,Y,9rsty\f,f)f[f° 



(2) 



where g rs is the metric tensor arising through the useage of 
curvilinear coordinates (3 03 ■ I n Cartesian coordinates g rs 
is presented by the unit matrix; T[, and XJ are two Jacobi 
marices. 

The tensor (|2jl is a quantitative measure of the deformation 
at given point of curvilinear coordinate system. Should a con- 
tinuous medium be represented as a collection of infinetise- 
mal cubes, then the tensor G^ completely describes how the 
edges and angles of these cubes contract, elongate and distort 
(as presented in Figure^. Differentiating the formula ©we 
arrive at the following evolution differential equation for the 
tensor Ga: 



dG 
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In order to describe plastic media, we separate the general 
deformation tensor into two tensors: G, and G are two parts 
of plastic deformation tensor of; Gk q is the tensor of elastic 
deformation: 



Gij — E E Gkq G? . 



(4) 



fc=l 9=1 



The choice of representing the plastic deformation tensor in 
two parts (instead of one as for elastic tensor) is related to the 
symmetry of the general tensor. The separation is the major 
step for the following consideration. The plastic deformation 
occurs in plastic materials that include pitch-like dense sticky 
liquids, solid amorphous materials, e. g. glasses, and poly- 
mers. These materials resist to the plastic deformation like 
solids and can flow like liquids. Unlike metals, where plastic 
deformation produces dislocations resulting in disorder of the 
crystalline grid, in our consideration, the plastic deformation 
does not change the material stucture. 

The separation produces separate evolution equations for its 
plastic and elastic parts. The evolution equation of the plastic 
deformation G can be presented as: 



dG 



k 3 



g L + J2v r VrG?=J2 d ^rV k 
r—1 r—1 
3 3 



(5) 



Here, as in equation Q v r is the velocity of medium; V is the 
spatial derivative; 9 is the physical parameter of continuous 
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medium associated with the relaxation rate of elastic defor- 
mation into plastic. 

Although the tensor Of is not symmetric, it can be obtained 
from the symmetric tensor 6ij through the standard procedure 
of index raising. Using equations Q, we derive the 
evolution equation of elastic deformation. The evolution of 
elastic deformation is an important factor in determining the 
stress of a medium. 



The parameter 9 has the units of inverse time: (sec" 1 ), and 
determines the stress relaxation rate. The external forces pro- 
duce stress and elastic deformation as a responce to this stress. 
In solely elastic material the responce is constant for constant 
deformation, and the parameter 9 is zero. In plastic materials 
9 is non-zero and the exponential reduction of elastic response 
is observed. For constant general deformation G, the elastic 
component G in (|4j decreases, while the plastic part G in- 
creases. This scenario is presented in Figure 2. 



^ + £V V r G kq = - J2 V k v r G rq 



J2 Gkr v,« r + k G rq + J2 G ^ K- 



(6) 



Equations Q and must be complemented by three bal- 
ance equations (three conservation laws for the flows of mass, 
momentum , and energy). The first is the mass balance equa- 
tion: 
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FIG. 2: Conceptual presentation of parameter 9. Elastic deformation 
produces stress relaxation resulting in plastic deformation through 
the conservation laws. 
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at 



(7) 
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where p is the density; v k is the velocity tensor. The next 
two equations are the momentum, and the energy balance 
equations: 



d(pv l 

at 



fc=i 



dfpH 2 \ ^„ k 



(8) 
(9) 



k=l 



Here /' is the density of external forces acting on the 
medium, e is the power of the forces; v , v 2 , v 3, are the com- 
ponents of velocity vector v; e is the specific inner thermal 
energy; H lk and w k are the momentum, and energy flows de- 
termined by the following formulae: 



Further analysis of this model indicates that for materials 
in which the plastic deformation does not change the struc- 
ture and produces no stress, it can excluded from consider- 
ation. This can be formulated through the '"forgetting prin- 
ciple'" . Suppose the medium evolves from an initial state at 
time t = 0, to the intermediate state, at t = to, and then 
continues evolution for t > to. The principle asserts that if 
the intermediate state total deformation of medium was solely 
plastic, then the medium evolution occurs without memory of 
the intermediate state. 

Equations 10-© form a complete set describing the evolu- 
tion of plastic medium in the frames of this model. They are 
used to describe the thermodynamics of plastic deformations. 
The specific internal thermal energy of the material depends 
on the entropy and on the elastic deformation. The differential 
of the energy can be written as: 
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where a lk is the stress tensor; 77 J9 is the viscosity tensor; 
x lk is the heat conductivity tensor and vj q is the tensor of 
velocity gradients: 



de = Td s -Y J Y j ~ at3dGl] 



i=i 3=1 



2p 



(10) 



Here, is an auxiliary tensor related to the stress tensor 
through the following formula: 



(ID 



Using equations dlOI . and (II Q the time derivative of e and 
its gradient can be found. Substituting them into equation (|9jl 
and using the equations (0-©, we obtain the equation for 
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specific entropy fl7il : 



Of 



k=l 



i=l 
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a%3 ® ij 
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(12) 



EEEE- 

i=l fe=l j=l g=l 



The three terms in the right hand side of equation d!2t 
describe three different mechanisms of entropy production: 
plasticity, heat transfer, and viscosity. Second and third terms 
of dl2l are typical for the visco-elastic solids. The term 



EE- 



i=i j=i 



T 
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responsible for plasticity is new one. It is non-zero indicat- 
ing the growth of entropy on the way to equilibrium. Elas- 
tic deformation produces stress leading to plastic deforma- 
tion, which in turn, results in additional entropy production 
as shown in d!2i . 

Once introduced, the material parameter 8 should be mea- 
sured experimentally. The parameter can be found using 
calorimetric measuremens, as can be seen from the tempera- 
ture dependence of the term dl 3i . Alternatively, this parameter 
can be measured through the static bending of a polymer, or 
dynamic flow of the material. These two possible approaches 
for 8 measurements are presented in Figure 3 below. 

In summary, we have suggested an approach for the sep- 
aration of plastic deformation within general nonlinear ten- 
sor deformation. This separation allows for a robust descrip- 
tion of the material's plasticity using the relaxation rate in- 
troduced through the parameter 8. This approach can be also 
used to describe experimental data of polymer deformation on 
nanoscale. 
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